Electronic correlations in transport through coupled quantum dots 
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The conductance through two quantum dots in series is studied using general qualitative ar- 
guments and quantitative slave-boson mean-field theory. It is demonstrated that measurements of 
the conductance can explore the phase diagram of the two-impurity Anderson model. Competition 
between the Kondo effect and the inter-dot magnetic exchange leads to a two-plateau structure in 
the conductance as a function of gate voltage and a two or three peak structure in the conductance 
vs. inter-dot tunneling. 
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The recent observation of the Kondo effect in trans- 
port through a quantum dot opened a new path for 
the investigations of strongly correlated electrons. Hav- 
ing confirmed earlier theoretical predictions j|]5), that a 
quantum dot behaves as a magnetic impurity, these ex- 
periments also serve as a critical quantitative test for ex- 
isting theories. In particular, unlike magnetic impurities 
in metals which have physical properties determined by 
the host metal and the impurity atom, the correspond- 
ing parameters in the quantum dot case can be varied 
continuously, enabling, for example, a crossover from the 
Kondo to the mixed-valence and the empty dot regimes 
in the same sample 

The behavior of a lattice of magnetic impurities, such 
as a heavy-fermion system, is characterized by the com- 
petition between the Kondo effect and the magnetic cor- 
relations between the impurities. An important step to- 
wards the understanding of this problem was taken by 
Jones and collaborators § , who studied the two-impurity 
problem. Their work demonstrated that this competition 
leads to a second-order phase transition when particle- 
hole symmetry applies. When this symmetry is broken, 
this transition is replaced by a crossover |t]-|9| . In view of 
the extensive experimental research on transport through 
two dots in series |l^Jll|], it is thus natural to try and 
understand how this phase-transition is manifested in 
the double-dot system, both because such systems may 
have important applications (such as a quantum-dot laser 
jUI), and because such a tunable system may reveal de- 
tailed information on the corresponding phase diagram. 

Transport through a double-dot system (see inset in 
Fig. has already received much theoretical attention, 
in particular in the high temperature, Coulomb block- 
aded regime 
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In experiments the Coulomb charg- 
ing energy and the excitation energy are much larger 
than temperature. Accordingly, only a single state on 
each dot is important, and double occupancy of each dot 
can be ignored. Denoting the energies of these states by 
(i = £o + and €2 = £o — Ae, respectively, and the 
tunneling amplitude between the dots by t, the isolated 



double dot system can contain zero, one or two electrons, 
depending on the chemical potential: N = for [i < e_, 
N = 1 for e_ < (i < e+, and N = 2 for \i > e+, with 
e± = eo ± \/Ae 2 + 1 2 . In the presence of a finite an- 
tiferromagnetic spin-exchange J between the dots, one 
still has the above three possibilities with e + replaced by 
e + — 3 J/4 |Q. In the Coulomb blockade regime there will 
be two peaks in the conductance vs. chemical potential 
at the degeneracy points /x = e±. Alternatively, starting 
from the N = 2 regime for t = 0, e+ will increase with 
increasing t, until it crosses the chemical potential and 
the ground-state will have a single electron in the double 
dot system. Again we expect a peak in the conductance 
at t — t + corresponding to £+(£+) = fx. 

At low temperatures, in addition to a renormalization 
of the above energy scales, e± — > e±, and t + — > t + , the 
Kondo effect starts to play a significant role in the trans- 
port. The relevance of the Kondo effect in the double 
dot system has been studied in 
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on the competition between the Kondo effect and anti- 
ferromagnetic exchange, and present detailed predictions 
for the conductance. (Recently Andrei et al. [|17[ have in- 
vestigated a very different realization of this competition, 
which applies to coupled metallic islands, close to points 
of charge degeneracy - i.e. near the Coulomb-blockade 
peaks |l8)). We find a rich phase diagram leading to in- 
teresting features in the conductance as a function of gate 
voltage and intra-dot tunneling. As the corresponding 
energy and temperature scales are experimentally acces- 
sible, these predictions are relevant to transport experi- 
ments in double-dot systems. 

To simplify notations, we assume in the following 
e\ = e 2 = eo and Vl = Vr = V, where Vl(Vr) is the 
coupling to the left (right) lead. Then the eigenstates of 
the double dot system are the even and odd states. As 
even-odd symmetry is broken anyway by the tunneling 
t, one can show that the above assumptions have little 
effect on the underlying physics. Following Ref. [^9|, the 
zero-temperature conductance per spin channel through 
the system can be expressed in terms of the retarded 
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Green functions for the even and odd states £/,T e *(aj) as: 

9 = ^T 2 \Q r e et {u = 0) - g r et {u = 0)| 2 where T = irpV 2 
and p is the density of states in the leads at the Fermi en- 
ergy. Defining the corresponding scattering phase-shifts, 
S a = 7r + argt/^ e * (/j,) (such that 5 a is in the 0-tt range), 
the conductance formula simplifies to 

g = — sin 2 (<5 e - S ) (1) 

The Friedel sum rule relates the total charge q per spin 
channel, on the double dot system, to these phase shifts: 
q = (S e + S )/tt. There is, however, no individual rela- 
tion between S a and the occupation of the corresponding 
state. 

For /i, > e_|_ (or alternatively t < EjT), there are N = 2 
electrons in the system (q= 1), and states with N = Q or 
N — 1 are high-energy states that can be eliminated from 
the Hilbert space. The effective low-energy Hamiltonian 
only involves spin degrees of freedom on the dots. It can 
be cast into the form of a model of two Kondo impurities, 
with Kondo couplings to the even and odd combinations 
of the conduction electrons in the leads, an inter- impurity 
magnetic exchange J oc t 2 /U , and a potential scattering 
term in the leads V e >° such that V e - V Q oc t V 2 /{el - 1 2 ). 

Previous studies of the two-impurity Kondo problem 
(mainly using Wilson numerical renormalization group 
(NRG)) J^,P have already yielded information on how 
the phase shifts S e , 5 a behave as a function of the cou- 
plings at T — in this regime. Let us start with the 
case J = 0. For t = there is an even-odd symmetry, 
and each channel has its own Kondo effect, leading to 
5 e = S = 7r/2, and naturally a zero conductance. As 
t/T is increased, the difference S e — S increases to reach 
the value 7r/4, at which point the conductance takes its 
maximum possible value: e 2 /h per spin channel. As t/T 
is further increased, the Kondo effect is gradually over- 
come by potential scattering and one reaches 5 ~ 0, 
S e ~ 7r in the large t/T limit (but still with t <C EjT). A 
slave-boson mean-field theory (SBMFT) presented below 
(cf. also Ref. jl6|]) yields in this regime: (J = 0, q ~ 1) 
S = S e — 5 a = 2 tan -1 t/T, leading to 

_ e 2 4t 2 r 2 

h {t 2 +T 2 ) 2 ( ' 

which reaches its maximum value at t = T (solid curve 
in Fig. The Kondo temperature in this regime is of 
order Tk = ci Tjj-e C2t / r , where the c's are weakly de- 
pendent on t/T and T% = iye _7r|£o|/r is the single-dot 
Kondo temperature. (The SBMFT yields ci = cos<S/2 
and C2 = (5/2) . The crucial content of that formula is 
that the coupled-dot Kondo temperature can be much 
larger than the single-dot Kondo temperature for small J 
and large t/T. This has important consequences for the 
observability of the effects described in this paper. 
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FIG. 1. Plot of the conductance vs. the tunneling between 
the dots, t, obtained by slave-boson mean-field theory. Due 
to the Kondo effect in the two-electron regime (N — 2) the 
conductance has a peak at t = V. As t increases beyond t+, 
the Kondo effect is quenched, until the ground-state contains 
a single electron (iV = 1), leading to a different Kondo state 
and an enhanced conductance. For finite antiferromagnetic 
coupling J (finite a — F 2 /UT^), the conductance peak is 
pushed to smaller values of t and becomes narrower, as the 
singlet formation destroys the Kondo state. In addition a 
second maximum in the conductance in the N — 2 regime 
emerges. Inset: The double dot system discussed in this pa- 
per. 

Let us now consider the effect of a finite J. For t = 0. 
the effective two-impurity Kondo model has particle-hole 
symmetry, and it is known from the work of Q that a 
phase transition exists at a critical value of the coupling 
J c /Tj£ ~ 2.2. For J < J c , the spin of each dot undergoes 
a Kondo effect with the leads and 5 e = 6 = ir/2. For 
J > J c , the two spins are locked into a singlet state and 
the Kondo effect does not apply, yielding 6 a — 0, 6 e = it. 
The phase-shift difference S jumps discontinuously from 
S = for J < J c to S = 7r for J > J c . (The conductance 
is, of course, zero for all J since t — 0). Turning on a 
small value of t/T is known to be a "relevant perturba- 
tion" on this critical point (with dimension 1/2, identical 
to that of J— J c ) |j],H and therefore smears the transition 
into a rapid crossover from 6 = to 5 = ir. For J close 
to J c , this smearing is described by a crossover scaling 
function: 

with (j)(x — > — oo) = and 4>(x — > +oo) = 1. As a re- 
sult, the conductance has a very sharp maximum as t/T 
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is increased for a fixed value of J close to J c . For J sig- 
nificantly larger than J c , the conductance remains very 
small with only a shallow maximum as t/T is increased. 
For intermediate values of t/T and J/T^, a quantitative 
calculation of 5 is needed in order to obtain the con- 
ductance, using e.g. NRG [§,|]2C]] or SBMFT. However, 
much can be said on a semi-quantitative level by using 
existing knowledge on the two-impurity Kondo problem. 
The phase shift 6 is an increasing function of J, which 
starts at the value given above Eq.(^), and increases un- 
til it saturates at 5 — tt at a scale J*. From the above 
estimate of Tjc the ratio J* /T^ increases exponentially 
with t/T. These considerations, and the knowledge of the 
crossover around J c (Eq.([|)), lead to a qualitative con- 
tour plot of the conductance in the (J/T^,t/T) parame- 
ter space, throughout the N ~ 2 regime, as displayed in 
Fig. | 




FIG. 2. Schematic contour map of the conductance in the 
N ~ 2 regime. Thicker lines denote higher conductance, the 
thickest one corresponding to g = e 2 fh. The broken lines are 
three physical contours (for different values of a = T 2 /UT K ) 
along which J ~ t 2 /U. 

In practice, the exchange J is not an independent 
parameter, but is a function of the interdot tunneling, 
J ~ t 2 /U . The contour plot above must thus be inter- 
sected by a curve J/T K = a(t/T) 2 , with a = T 2 /UT^, 
in order to follow the dependence of the conductance as 
a function of t/T. Since is a very sensitive function 
of the energy scales (such as eo and T), the control pa- 
rameter a can be varied continuously over many orders 
of magnitude, allowing an experimental investigation of 
most of the phase-diagram. Thus, as a function of t, the 
maximum conductance e 2 /h is reached for t ~ T with a 
peak width At oc T for small a, while the peak is pushed 
down to much lower transmission t ~ T/^/a and becomes 
very narrow At ~ T/a for large a. In addition, as the 



saturation scale J* increases exponentially with t, one 
may expect, for an intermediate a (middle broken curve 
in Fig. ||), an additional peak in the conductance vs. t 
in the N — 2 regime. These results are indeed confirmed 
by the SBMFT calculation (see Fig. [[J). 

As t is further increased (t > the equilibrium 

charge decreases to N = 1 (q = 1/2). In this regime the 
effective Hamiltonian is that of a single-impurity Kondo 
problem in the even parity sector Jl5| , leading to unitary 
scattering S e ~ 7r/2. In the odd parity sector, we have 
an almost empty resonant level with 5 a ~ (Note that 
(<5 e + S )/tt = q ~ 1/2 consistent with Friedel sum rule). 
Throughout this regime, we therefore expect the zero- 
temperature conductance to be maximum g = e 2 /h and 
essentially independent of t. In this regime, the inter-dot 
exchange J plays little role. 

Similar interesting behavior is expected as a function 
of gate voltage, that controls the depth of the level energy 
eo with respect to the chemical potential (see Fig. ||) . For 
a very deep level the Kondo temperature is exponentially 
small, and thus J/T^ is large and quenches the Kondo 
effect. As eo increases, the Kondo temperature increases 
and we enter the (N = 2) Kondo state, and a finite con- 
ductance. This conductance remains constant (at zero 
temperature) at a value smaller than e 2 /h, determined 
by the value of t, until eo crosses the Fermi energy and 
a new (N — 1) Kondo state is formed. There the con- 
ductance is given by its maximum value, e 2 /h per spin. 
As eo is further increased the double dot system becomes 
empty and the conductance drops to zero. 
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FIG. 3. Plot of the conductance vs. the level energy, as 
obtained from SBMFT (for t = 2 and U/T = 10 4 ). The con- 
ductance rises horn a very small value (the singlet regime, 
J 3> Tk), to a independent value (TV = 2 Kondo regime, 
J <JC Tk), and then to g = e 2 /h (N — 1 Kondo regime) 
before dropping to zero for an empty dot. 
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To substantiate these semi-quantitative arguments we 
have performed a quantitative calculation of the phase 
shifts and the conductance using a slave-boson mean-field 
approximation. This method becomes exact as the num- 
ber of spin-degrees of freedom goes to infinity, and has 
been previously used in order to study the two-impurity 
Anderson model in Ref. raj. It was recently applied in the 
present context in Ref. U& but only in the case J = 0. 
We have solved numerically the full set of equations in- 
cluding the tunneling t and exchange J, but we only 
quote here the simplified version of the equations |7j that 
hold in the q = 1 regime (one electron in the double dot 
system per spin-state, corresponding to N — 2). For 
small enough values of J/T^, the phase-shift difference 
S is given by the solution of: 



2tt 



M/2T 



t s 

— COS 

I 2 



J 



(4) 



As J is increased beyond a critical coupling Jf s , S 
reaches the value it: this is either a smooth transition 
for t > 1/tt, or a first-order jump for t < 1/n (deter- 
mined by free-energy considerations). The existence of 
a phase transition even for non-zero values of t/T is an 
artifact of the SBMFT approximation: J^ B should ac- 
tually be interpreted as an estimate of the saturation 
scale J* discussed above. This spurious transition does 
not affect qualitatively the behavior of the conductance, 
except when it becomes very small: there a strictly zero- 
value of g can be found (as evident on Fig. [l]) whereas 
the real system would have only a very small but finite 
g. The SBMFT also provides a quantitative estimate of 
the Kondo scale for the coupled dot system in the q ~ 1 
regime, as mentioned after Eq. (||). 

In conclusion, we have demonstrated that measure- 
ments of the conductance through a double-dot system 
can explore the phase diagram of the two-impurity An- 
derson model. By changing the control parameter a — 
r 2 /UTft (which depends sensitively on gate- voltage) , 
one can make various cuts through the phase-diagram 
(Fig. |^), leading to non-trivial features in the conduc- 
tance vs. gate-voltage and inter-dot tunneling (Figs, [j] 
and ^]). As the relevant temperature scale can be much 
higher than the single-dot Kondo temperature we believe 
that these predictions could be tested experimentally. 
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